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Studies  on  the  Gas  Bubble  Resulting  From 
Underwater  Explosions 

ON  THE  BEST  LOCATION  OP  A  MIKE  NEAR  THE  SEA  BED 

At  the  request  of  the  David  Taylor  Model  Basin  of 
the  Bureau  of  Ships  and  of  Division  Re  2  of  the  Btireau  of 
Ordnance,  a  series  of  investigations  has  been  undertaken 
by  the  New  York  University  Group  of  the  Applied  Mathematics 
Panel  of  the  N.D.R.C.  with  the  object  of  analyzing  the 
phenomena  associated  with  the  gas  bubble  produced  by  an 
underwater  explosion* 

The  present  study,  carried  out  by  Dr.  Max  Shiffman 
and  Dr.  Bernard  Friedman,  with  the  cooperation  of  the 
Mathematical  Tables  Project  in  the  extensive  numerical  work, 
is  concerned  with  the  following  problem:   If  a  mine  of  given 
weight  of  explosive  is  to  be  placed  near  the  sea  bed,  what 
position  should  it  have  to  cause  maximum  damage  to  a  target 
at  the  surface? 

It  seems  appropriate  to  make  a  few  introductory  re- 
marks about  the  broader  research  program  of  which  this  report 
is  a  part.  Often  the  destructive  effect  of  an  underwater 
explosion  is  not  wholly  due  to  the  high  pressure  shock  of  the 
explosion.  After  the  initial  shock  wave  has  passed  with 
enormous  speed,  a  comparatively  slow  pulsation  of  the  gas 
bubble  (consisting  of  the  biirned  gases)  takes  place.   In  the 
second  and  sometimes  even  third  pulse  of  this  motion  of  the 
bubble,  a  pressure  pulse  of  considerable  strength  is  emitted. 
While  the  peak  of  the  shock  pressure  from  the  explosion  is 
far  above  tlxat  of  the  later  pulses  (in  typical  cases,  six 
to  ten  times  the  pressure  of  the  second  pulse),  the  duration 
of  these  later  pulses  is  much  longer  (about  twenty  times  as 
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long);  so  that  the  moraentujii  imparted  and  the  destruction 
caused  may  be  comparable  to  that  of  the  first  shock.   The 
importance  of  the  second  pulse  has  been  confirmee'  more  and 
more  not  only  by  model  experiments,  but  also  by  analysis 
of  actual  damage  to  allied  ships  by  mines  [1],  [2]« 

Just  as  significant  as  the  long  duration  of  the  sec- 
ond press\ire  pulse,  is  the  migration  of  the  gas  bubble 
during  the  pulsation.   The  migration  can  increase  the  effect 
of  the  second  pressure  pulse  because,  under  suitable  conditions, 
the  center  of  this  pulse  might  be  much  closer  to  the  target 
than  the  explosion.  Factors  influencing  the  motion  of  the 
bubble  are:  rigid  walls  such  as  the  hull  of  a  ship  or  the 
sea  bed,  which  attract  the  bubble;  the  free  surface  of  the 
water,  which  has  a  repulsive  effect;  and  the  buoyant  force 
of  gravity,  which  causes  the  bubble  to  rise  tov/ards  the  sur- 
face. 

Understanding  and  controlling  the  complicated  inter- 
play of  these  effects  is  of  sufficient  practical  importance 
to  warrant  comprehensive  research.   Unfortunately,  a  great 
variety  of  experiments  with  full  size  charges  is  hardly 
feasible,  and  experiments  with  small  charges  cannot  easily 
be  reinterpreted  for  large  charges.   Such  difficulties  call 
for  theoretical  investigation  as  a  guide  to  experiment. 

Historically,  the  earliest  theoretical  problem,  satis- 
factorily analyzed  in  older  studies,  is  that  of  a  single 
spherical  gas  bubble  immersed  in  an  infinite  body  of  water; 
in  particular,  the  pressure  as  a  function  of  time  and  distance, 
and  the  period  of  pulsation  of  the  bubble  can  be  explicitly 
determined. 

For  the  past  few  years  an  extensive  research  program 
on  underwater  explosions  has  been  pvirsued  by  G.  I.  Taylor 


•K-  Bracket  references  refer  to  the  bibliography. 
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and  his  associates  in  England.   In  this  country,  theoretical 
work  was  started  by  C.  Herring  (Division  6,  N.D.R.C  )  and 
is  being  carried  on  at  the  David  Taylor  Model  Basin  (in  par- 
ticular by  Captain  W.  P.  Roop,  E.  H.  Kennard,  G.  E.  Hudson) 
and  by  the  N.Y.TJ. Group  of  the  A.M. P.  at  the  request  of  the 
Model  Basin  and  of  the  Bureau  of  Ordnance;  considerable  experi- 
mental research  is  being  done  at  the  Underwater  Sound  Labora- 
tory of  Division  8,  N.D.R.C 

Taylor  [3]  and  others  have  determined  the  upward  motion 
of  the  bubble  under  the  buoyant  force  of  gravity,  a  motion 
taking  place  in  jerks.  Herring  [4]  and  others  have  studied 
the  manner  in  which  rigid  walls  attract  the  bubble,  assuming 
that  the  bubble  is  rather  far  from  the  wall.   Shiffraan  [5] 
has  developed  an  improved  method  which  permits  numerical 
analysis  of  the  motion  when  the  bubble  is  close  to  the  wall, 
and  even  when  it  touches  the  wall.   The  influence  of  non- 
rigid  walls  and  the  repulsive  effect  of  the  free  surface  of 
the  water  on  the  bubble  have  also  been  studied.   Much  of  the 
material  is  condensed  and  supplemented  in  the  comprehensive 
work  by  Kennard[6]. 

A  very  important  aspect  of  the  problem  is  concerned 
with  plastic-elastic  deformations  of  the  target  in  inter- 
action with  the  motion  of  the  water.   The  process  of  damage 
can  be  understood  only  by  a  careful  analysis  of  this  inter- 
action.  Extended  research  by  J.  G .  Kirkwood  [7]  has  estab- 
lished that  the  damage  to  a  structure  depends  on  the  ratio  of 
the  duration  time  of  the  incoming  pressure  pulse  to  the  "time 
constant"  of  the  structure,  i.e.,  the  span  of  time  during 
which  the  structure  is  "receptive"  to  the  impinging  pressure 
effects.   If  this  ratio  is  small,  the  impulse  is  the  more 
Important  factor  for  damage,  while  if  it  is  large,  the  peak 
pressure  is  more  important. 

Another  problem  that  should  be  mentioned  is  that  of 
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the  change  In  the  shape  of  the  gas  bubble.   The  migrating  gas 
bubble  does  not  retain  its  spherical  shape,  but  flattens  out 
in  a  plane  perpendicular  to  the  direction  of  the  motion. 
The  result  is  a  deceleration  of  the  bubble  and,  as  a  conse- 
quence, an  intensification  of  the  pressure  pulse.   Prelim- 
inary studies  of  this  effect  have  been  carried  out  independently 
in  England  and  at  N.Y.U.  by  Dr.  Shiffman  [8]. 

Finally,  the  still  unsolved  problems  presented  by  the 
striking  phenomena  occurring  at  the  water  surface  must  be 
mentioned.  A  proper  interpretation  of  the  "domes"  and 
"pliimes"  rising  into  the  air  above  the  explosion  should  yield 
much  information  concerning  the  process  taking  place  under 
water. 

The  effect  of  gravity  in  moving  the  bubble  closer  to 
the  target  is  not  the  only  factor  which  leads  to  increased 
damage  by  the  second  pulse.   Another  element  of  major  im- 
portance, emphasized  by  Shlffman  and  Friedman,  is  that  the 
peak  pressure  of  the  second  pulse  depends  greatly  on  the 
state  of  motion  of  the  bubble  at  the  moment  of  the  first 
contraction.   It  is  shown  that  the  peak  pressure  of  the  sec- 
ond pulse  possesses  a  decided  maximun  if  the  gas  bubble  is 
stationary  at  the  moment  of  contraction.   For  depths  of 
water  and  weights  of  explosives  actually  used,  it  is  indeed 
possible  to  keep  the  bubble  stationary  by  counterbalancing 
the  gravitational  force  with  the  attractive  force  of  the 
sea  bed  and  the  repulsive  force  of  the  water  surface. 

This  principle  of  stabilization  forms  the  main  point 
of  the  following  report.   It  is  assumed  that  the  mine  is  to 
be  placed  fairly  close  to  the  sea  bed,  but  not  necessarily 
directly  on  it,  and  it  is  then  found  (roughly  stated)  that 
for  mines  of  about  500  to  2000  lbs.  of  explosive  in  water 
about  70  to  150  ft.  deep,  the  optimal  position  is  about  one 
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maximiun  bubble  radius  above  the  sea  bed;  this  means  tliat 
the  mines  should  be  moored  approximately  20  ft.  above  the 
sea  bedo   More  precisely,  for  water  of  given  depth  and  for  a 
mine  of  given  weight  of  explosive,  there  is  a  definite  opti- 
mum location  for  the  detonation,  that  location  being  given 
by  the  principle  of  stabilization. 

Of  course,  such  a  statement  makes  sense  only  if,  at 
the  outset,  the  position  of  the  mine  Is  restricted  to  a 
place  deep  under  water,  i.e.,  fairly  distant  from  the  target. 

Similar  problems  concerning  the  optimal  position  of 
an  explosion  also  occur  for  charges  nearer  the  surface,  and 
for  charges  dropped  from  fast  flying  planes  and  thus  having 
an  initial  horizontal  and  vertical  momentu:n  at  the  instant 
of  the  detonation.   The  discussion  of  such  questions  is 
planned  for  later  reports. 

Richard  Courant 
Contractor's  Technical  Representative 
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Part  1«  Discussion  of  the  Restilts 

lo  Preliminary  Remarks « 

The  present  report  is  concerned  with  the  following 
problem:   what  Is  the  "best  location  of  a  mine  near  the  bot- 
tom of  the  sea  so  that  dujrlng  the  secondary  pulse  the  peak 
pressure  Is  a  maximum?   The  result  obtained  Is:  for  weights 
of  500  to  2000  lbs.  of  T.N.T. ,  In  water  of  depth  100  to 
150  ft. ,  the  mine  should  be  placed  about  20  ft.  above  the 
sea  bed.   More  exact  locations  are  contained  In  table  1, 
p.  9. 

It  is  significant  to  compare  different  mine  positions 
as  to  their  effects  on  the  secondary  pulse.   This  report  demon- 
strates that,  during  the  secondary  pulse,  the  highest  possible 
peak  pressure  at  the  surface  Is  obtained  when  the  mine  Is 
placed  in  a  definite  "best"  position.   The  theor'y  indicates 
that  if  the  mine  is  placed  even  a  few  feet  away  from  its 
best  location,  the  peak  pressure  might  be  as  little  as  one- 
half  of  that  corresponding  to  the  "best"  location."  Further- 
more, when  the  mine  is  placed  directly  on  the  sea  bed,  experi- 
mental evidence  is  available  shov/lng  that  the  secondary  pulse 
is  weak  and  erratic. [9] 

Another  point  of  interest  is  that  the  "best"  position 
also  Increases  the  effectiveness  of  the  shock  wave  compared 
with  a  location  on  the  sea  bed  itself.   In  the  latter  case, 
some  experiments  Indicate  that  the  presence  of  the  sea  bed 
has  no  effect  on  the  shock  wave,  and  others,  that  it  results 
in  an  Increase  of  about  12  percent  in  the  peak  pressure  and 

%   See  the  example  discussed  in  section  4  below. 
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impulsed   On  the  other  hand,  when  a  mine  is  at  the  "best" 
location  (about  20  ft.  above  the  sea  bed  for  a  depth,  say, 
of  100  fto),  then  the  shock  wave  peak  pressure  at  the  sur- 
face is  increased  about  25  percent  over  that  of  the  sea  bed 
position  since  the  pressure  at  the  svirface  is  inversely 
proportional  to  the  distance  from  the  explosion. 

2«  The  principle  of  stabilization. 

The  best  location  of  the  mine  is  determined  by  the 

following  Principle  of  Stabilization:  For  a  given  mass  of 

e?:plo3ive,  the  maxinxun  peak  pressure  in  the  secondary  pulse 

is  obtained  if  the  gas  bubble  produced  by  the  explosion  is 

kept  motionless  at  the  time  of  its  minimum  size.   This 
stabilization  can  be  attained  for  mines  by  suitably  balancing 

the  upward  buoyant  force  on  the  bubble  against  the  dov/nward 

attractive  force  due  to  the  sea  bed.   In  other  cases,  the 

repulsive  force  of  the  surface,  as  well  as  additional  factors, 

could  be  utilized. 

In  part  II,  a  mathematical  demonstration  of  the  principle 

of  stabilization  is  given.   The  following  plausible  argument 

also  points  to  the  same  result.   In  general,  the  total  energy 

of  the  bubble  and  water  is  divided  into  two  parts,  namely, 

the  kinetic  energy  of  the  water  surrounding  the  bubble,  and 

the  internal  energy  of  the  gas  inside  the  bubble.   If  the 

bubble  is  moving  at  the  time  of  its  minimum  size,  then  some 

of  the  total  energy  is  diverted  into  kinetic  energy  of  the 

water  and  thus  less  is  available  for  the  internal  energy 

Vc   [9]  reports  no  change  in  the  shock  wave  for  a  mud  bottom. 
Experiments  in  Woods  Hole  for  a  sand  bottom  show  an  increase 
of  12  percent,  according  to  an  oral  discussion.   Theoretically, 
if  the  sea  bed  were  perfectly  rigid,  the  peak  pressure  would 
be  multiplied  by  the  factor  -^g"  ,  which  means  an  increase  of 
26  percent.   Actually,  however,  the  explosion  tears  a  hole  in 
the  sea  bed  and  also  transmits  a  shock  into  ito 
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of  the  gas.  This  would  reduce  the  pressure  inside  the  gas 
and  also  the  pressure  in  the  water.  To  obtain  the  optlraun 
pressiiTe,  therefore,  this  bubble  must  be  motionless  at  the 
time  of  its  miniimira  size." 

For  an  appraisal  of  the  value  of  stabilization,  the 
damage  due  to  the  secondary  pulse  must  be  studied.   All 
theories  of  damage  show  that  for  a  pressure  pulse  which 
lasts  long  relative  to  the  time  constant  of  the  structiire 
to  be  damaged,  the  damage  is  approximately  proportional  to 
the  peak  pressure  and  not  to  the  impulse  or  energy.   Since 
the  duration  of  the  secondary  pulse  is  approximately  twenty 
times  that  of  the  shock  wave,  it  would  seem  that  maximizing 
the  peak  pressure  of  the  secondary  pulse  increases  the  damage. 

3.  The  location  of  the  mine. 

Let  'vY  be  the  weight  of  the  explosive  in  pounds, 
(the  numerical  values  refer  to  T.N.T.  for  which  the  charac- 
teristic constants  were  available  to  the  authors),  D  the 
distance  in  feet  from  the  bottom  of  the  sea  to  a  point  33 
feet  above  sea  level  (allowing  for  the  pressure  of  the  atmos- 
phere), and  B   the  distance  in  feet  from  the  center  of  the 
mine  to  the  sea  bed.   The  problem  is  to  determine  the  best 
value  for  B   if  W  and  D  are  given.   This  is  more  easily 
expressed  in  terms  of  the  non-dimensional  quantities  which 
are  used  in  part  II.   For  the  unit  of  length  select 

/  y^ 
(1.1)  L  =  13.2  gj   ft., 

which  represents  (approximately)  the  maximum  radius  of  the 

■K-  This  argument  is  not  completely  rigorous  since  contri- 
butions to  the  pressure  in  the  water  occur  from  the  motion 
of  the  water  as  well  as  from  the  gas  pressure  inside  the 
bubble.   The  fact  that  the  gas  pressure  is  the  more  important 
contribution  arises  from  the  mathematical  analysis  of  the 
interplay  of  the  tv/o  effects. 

CONFIDENTIAL 


CO^IFIDENTIAL 


-9- 


gas  bubble  at  the  first  expansion. *  Set 


(1.2) 


H  -  D   >,  _  ^o 


so  that  d  and  b   represent  the  distance  D  and  B 
measured  in  "bubble  radii". 


t 
33  ft. 


D  ft 


T 


Sea.    leueL 


'  ^w  '////////  ^  /  ^/  y  ^  ///  / 


The  best  position  of  the  mine  is  determined  approxi- 
mately from  the  equation 


(1.3) 


d  =  6.2b^  +  3.3b  +  .4 
o       o 


(See  derivation  in  part  III).   Graphs  of  (1.1)  and  (1.3) 
are  drawn  in  figures  1,  2  at  the  end. 

On  the  basis  of  (1.1)  and  (1.3),  the  following  table 
can  be  constructed  showing  the  best  location: 

Table  1 


Weight  of  T.N.T. 
(in  lbs. ) 

Depth  of  Sea  Bed 

From  Sea  Level 

(in  ft.) 

Best  Location  of  Mine 
from  Sea  Bed 
(in  ft.) 

1500 

150 

21.0 

1000 

150 

20.2 

1000 

100 

17.0 

500 

100 

15.5 
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The  following  limitations  and  qualifications  should 
be  noted.   First,  the  effect  of  the  free  s\irface  has  been 
ignored,  so  that  these  results  are  valid  if  the  mine  is 
sufficiently  far  from  the  sea  level,  say,  more  than  3.5  L  ft. 
Second,  it  is  assumed  that  the  bubble  does  not  meet  the 
sea  bed  in  the  couj?se  of  its  motion.   This  is  approximately 
equivalent  to  B  >  .SL.   Third,  the  sea  bed  is  assumed 
to  act  like  a  rigid  wall. 

The  validity  of  this  last  assizmption  should  be  tested 
by  experiment.  But  there  are  plausible  reasons  for  supposing 
that  as  far  as  the  balancing  phenomenon  is  concerned,  the 
sea  bed  does  act  very  much  like  a  rigid  wall.   During  the 
largest  part  of  the  time  of  pulsation  of  the  bubble,  the 
bubble  is  large  and  the  pressure  in  the  water  is  low;  there- 
fore, a  sand  bottom  could  be  considered  as  rigid.   The 
balancing  effect  is  determined  in  the  main  by  thiK  portion 
of  the  period  of  pulsation. 

4.  An  exaii^le. 

Consider  a  mine  with  1500  lbs.  of  T.N.T,  in  water 
of  depth  150  ft.   Figure  3  is  a  graph  showing  the  relative 
magnitudes  of  the  peak  pressujre  in  the  secondary  pulse 
when  the  mine  is  placed  at  varying  distances  from  the 
bottom  of  the  sea.   The  best  location  of  the  mine  is  21.0  ft. 
from  the  bottom;  the  resulting  peak  pressure  at  any  point 
in  the  water,  as  computed  in  part  II,  is  2015/R  atmospheres, 
where  R   is  the  distance  in  feet  from  the  point  to  the  center 
of  the  contracted  bubble.   If  the  mine  is  initially  placed 
at  any  other  distance  from  the  bottom,  this  value  for  the 
peak  pressure  is  to  be  multiplied  by  the  factor  whose  graph 
is  drawn  in  figure  3. 

This  graph  demonstrates  a  remarkable  sharpne s s^"* in  the 
peak  pressure  curve  as  a  function  of  the  distance  of  the 

-::-  The  correction  due  to  the  free  siorface  is  treated  in 
part  III,  section    S^ 

-x-  In  reality  the  graph  will  not  be  quite  as  sharp  because 
of  deviations  from  the  assxxmptlons  listed  on  page  12. 
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mine from  the  bottom.   For  example,  if  the  mine  Is  placed 
either  17  ft.  or  29  ft.  from  the  bottom,  the  peak  pressure 
Is  only  one -ha  If  of  the  maximum. 

If  one  is  interested  in  the  damage  to  objects  on  the 
surface  of  the  water,  it  is  necessary  to  take  into  account 
the  migration  of  the  gas  bubble.   Again  considering  the  above 
exaiaple,  the  peak  pressure  at  the  surface  directly  overhead 
is  15.6  atmospheres  (excluding  the  reflection  from  the  bottom) 
when  the  mine  is  placed  21.0  ft.  from  the  bottom.  For  any 
other  location  of  the  mine,  this  value  is  to  be  multiplied 
by  the  factor  whose  graph  is  drawn  in  figure  4.   The  graph 
still  exhibits  a  considerable  sharpness,  but  if  one  places 
the  mine  somev/hat  higher  tlian  21.0  ft.  from  the  sea  bed,  the 
effect  on  the  surface  would  be  improved  slightly.   This  is 
due  to  the  upward  migration  of  the  bubble.   Note,  however, 
that  if  the  mine  is  placed  too  high,  say  27  ft.  from  the 
bottom,  the  favorable  effect  of  the  upward  migration  is 
sharply  counterbalanced  by  the  resulting  weakness  of  the 
secondary  pulse. 

Part  II.  Mathematical  Study  of  the  Secondary 
Pressure  Pulse 

1.  Introduction. 

High  pressure  pulses  are  produced  only  when  the  size 
of  the  gas  bubble  is  near  its  minimum.   During  this  time, 
the  buoyant  force  due  to  gravity  is  small.   Likewise,  the 
proximity  of  rigid  walls  or  free  surfaces  will  not  materially 
affect  the  motion,  since  their  influence  depends  on  the  ratio 
of  the  bubble  radius  to  the  distance  from  the  bubble,  and 
is  small  if  this  ratio  is  small.   Therefore,  during  the  stage 
of  minimum  size,  the  bubble  can  be  considered  as  immersed 
in  an  infinite  body  of  water  and  subject  to  no  outside  forces. 
The  pressure  pulse  produced  by  the  bubble  under  these 
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clrcuinstances  will  be  investigated. 

On  the  other  hand,  all  these  outside  influences 
do  affect  the  bubble  in  the  stage  when  it  is  not  small. 
Consequently,  in  contracting  to  its  minimum  size, the 
bubble  and  surrounding  water  have  a  linear  momentum  which 
remains  practically  constant  while  the  bubble  passes  through 
the  minimum  size  stage  when  the  outside  influences  can  be 
neglected.   The  pressure  pulse  produced  by  the  bubble  depends 
on  the  linear  momentum  acquired  by  the  bubble.   The  dependence 
of  the  peak  pressure  on  this  linear  momentum  is  the  main  objec- 
tive of  this  investigation  and  is  represented  in  figure   b- 

For  the  derivation  of  this  result,  the  following  assump- 
tions are  made: 

!•  The  water  is  an  ideal  incompressible  fluid. 

2»  The  bubble  remains  spherical  in  shape. 

3.  The  gas  inside  the  bubble  is  in  thermal 

equilibrium  at  each  instant  and  follows  the 
adiabatic  law. 
These  assumptions  are  reasonable  for  the  major  portion  of 
the  period  of  pulsation,  and  are  violated  only  in  the  very 
short  time  interval  when  the  radius  of  the  bubble  is  small. 
The  violations  are  in  the  nature  of  corrections  to  the  theory, 
and  will  not  materially  affect  the  principle  of  stabilization. 
We  shall  not  enter  into  a  discussion  of  the  assumptions,  but 
merely  refer  to  [4]. 

2.  The  Energy  Equation. 

Let  A  be  the  radius  of  the  bubble  at  any  time,   B 
the  vertical  distance  of  its  center  from  some  horizontal 
level,  and  P   the  hydrostatic  pressure  of  the  water  at  the 
center  of  the  bubble.   The  motion  of  the  bubble  is  described 
by  specifying  A,  B  as  functions  of  the  time   T» 
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/^/////  /////?//  ^  /  /  /  /  /    /// 


The  velocity  potential  $  describing  the  flow  around 
a  moving,  pulsating  sphere  in  an  infinite  body  of  water  is 
easily  obtained  frora  classical  hydrodynamics.   It  is 


(2.1) 


<!>  = 


_  A^A»    A^B'   cos  9 
R  ^  "2  ^2"  *' 


where  R,©  are  coordinates  as  indicated  in  the  diagram  above, 
and  the  "prime"  denotes  a  time  derivative.  The  kinetic  energy 
Y  of  the  water  is  then 


(2.2) 


f     =  2rf/)A^[A'^  +  iB'^] 


where  p  is  the  density  of  the  water.  For  these  classical 
results,  see,  for  example,  [3]  or  Appendix  II. 

The  potential  energy  XL    of  the  system  of  gas  bubble 
and  water  is 


(2.3) 


U    =   ^ITA^P^  +  G(A) 


where  the  first  term  is  the  gravitational  potential  energy 
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of  the  bubble  due  to  the  absence  of  water  and  the  second 
term  G(A)   is  the  Internal  energy  of  the  gas.   If  the 
gas  obeys  the  adlabatlc  law,  the  Internal  energy  is 

/  N     K  M^  ._ 
G(A)  =  -3(^-17 

A 

where  i      is  the  adlabatic  exponent,   K   is  a  constant  depend- 
ing on  the  explosive,  and  M  is  the  mass  of  the  explosive 
used. 

The  energy  equation  is 

(2.4)  r|^+U=E 

where  E  is  the  constant  total  energy  in  the  system  (after 
the  passage  of  the  initial  shock  wave  due  to  the  explosion). 

3»  Non-dimensional  variables. 

A  considerable  simplification  in  the  writing  of  the 
equations  is  obtained  by  introducing  non-dimensional  quantities 
with  appropriate  scaling  factors.   Likev/ise,  in  part  III,  it 
is  even  more  important  to  introduce  non-dimensional  quantities. 
For  the  purpose  of  coinparing  the  formulas  developed  in  this 
part  with  those  in  part  III,  we  shall  use  the  scaling  factors 
convenient  for  part  III.   It  should  be  mentioned  that  these  are 
not  the  most  convenient  factors  to  use  if  the  results  of  part 
II  were  the  sole  objective. 

Set 

(2.5)  A  =  La,   B  =  Lb,   T  =  Ct, 

where   L,C  are  scaling  factors  with  the  dimensions  of  length, 
time, respectively,  and  a,  b,  t  are  non-dimensional  variables. 
(Henceforth,  capital  letters  will  indicate  dimensional  quantities 
and  small  letters  non-dimensional  quantities.)   The  scaling  factors 
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which  we   shall  use  are 

(2.6) 

C   =  L 


3E 


where  P   Is  the  hydrostatic  pressxire  at  the  center  of 
the  explosion,  and  p   is  the  density  of  water*'*"  These 
scaling  factors  are  advantageous  because  they  have  a  simple 
physical  interpretation  intimately  connected  with  the  motion 
of  the  bubble.   It  will  be  seen  in  part  III  that  L  is  the 
maximum  radius  of  the  bubble  if  the  internal  energy  of  the 
gas  could  be  neglected.   Actually  the  maximum  radius  is 
approximately  •92L.   Also,  it  will  be  shown  in  part  III  that 
the  time  scaling  factor  C  represents  2/3  of  the  period  of 
pulsation  of  the  bubble  if  no  outside  influences  such  as  rigid 
walls  or  free  surfaces  are  present. 

The  energy  equation  in  terms  of  the  non-dimensional 
variables  becomes 


(2.7)  a2p2.1J2].a=.0^=X, 


where  the  "dot  "  denotes  differentiation  with  respect  to  the 
non-dimensional  time   t. 

4,  The  scaling  factors  and  the  internal  energy  of  the  gas  for 

T.  N»  T. 

It  is  desirable  to  express  the  scaling  factors   L,C 
in  terms  of  the  mass  of  explosj/^e  and  the  depth  of  the  explosion. 


«•  These  scaling  factors  differ  from  those  used  by  G.  I.  Taylor 
in  [  3] . 
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Let  D   be  the  distance  of  the  center  of  the  explosion  from 
o 

a  point  33  ft.  above  sea  level,  so  that  P   =  pgD  •   It  is 
still  necessary  to  express  the  energy  constant  E  in  terms 
of  the  mass  of  explosive*   For  the  case  of  T,  N.  T» ,  experi- 
mental results  indicate  that  E   is  approximately  one-half  of 
the  total  chemical  energy  released  by  the  explosion.   See  [4]. 

Slight  changes  in  E  will  not  materially  affect  the  results 

3 1— 
since  E   appears  in  the  form   \/e^   Using  the  values  given 

by  G.  I.  Taylor  in  [  3] ,  p.  4 ,  and  converting  to  the  English 

system  of  measurements,  one  finds  that 

L  =  13.2  W^  feet. 


(2.S) 


wV3 
C  =  2»85  — ^^   seconds  , 

o 


where  W  is  the  weight  of  T.  N.  T.  in  pounds,  and  D  is 
the  distance  in  feet  of  the  center  of  the  explosion  from  a 
point  33  ft.  above  sea  level. 

Likewise,  making  use  of  the  experimental  values  given 
by  G.  I.  Taylor  in  [3],  p.  4  ,  the  quantity  ^\^^   for  T.  N.  T. 
is,  in  metric  units, 

0(A)    I'lSSM^y^ 
-E 737^   ' 


where     M  is    the  mass   of  T,    N,    T.    in  grams.      Changing   to 

gr 

non-dimensional  variables,    one   obtains 


(2.9) 


G(A)    _       k 


1.  183  M     ^^^ 

where     k     =  -^^     < 

L      3/4 
cm 
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Changlng  to  English  units  and  using  (2»8),  one  finds 
that 

(2.10)  k  =  .0607  dJ        , 

where  D   has  the  same  meaning  as  in  (2«8).   Equations 
(2«9)  and  (2.10)  are  the  desired  expressions  of  the  internal 
energy  for  T«N.  T. 

5«  The  equations  of  motion  for  small  values  of  the  radius  a» 

The  secondary  high  pressure  pulse  is  emitted  when  the 
bubble  is  near  its  minimum  size,  and  henceforth  we  shall 
consider  the  case  when  the  radius  is  small.   The  additive 
term  a   in  the  energy  equation  (2.7)  is  then  negligible 
compared  with  the  number  1,  and  the  energy  equation  becomes 

(2.11)  a3(a2  +  I  b^)  +  -^  =  1. 
The  Lagrangean  equations  of  motion  are 

<^-i^'  1  If  -  If  =  ° 

where  ^=  j  -  "U.  •   They  become 

(2.14)  1^  (2  a^a)  =  3a2(a2  .  i  b')  -  f  ;^   , 

(2.15)  |^(|  a^  b)  =  0. 
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The  equation  (2. 15)  shows  that  when  the  bubble  is  near 

1  3» 
its  minimum  size,  the  linear  momentum  t  ^  t>   remains  constant. 

Denote  this  constant  value  by  s  .   Then       ; 

(2.16)  b=.^   , 

and  the  energy  equation  (2.11)  yields 

-2 

The  equations  (2.16),  (2.17)  completely  determine,  in  terms 
of  s,  the  motion  of  the  bubble  during  its  most  interesting 
phase  when  the  radius   a   is  small.   An  integration  of  (2.17) 
and  then  of  (2.16),  yields   a,  b  as  functions  of   t. 

Equation  (2.14)  gives  an  expression  for  the  quantity 
(a  a)*,  which  is  important  for  the  determination  of  the  pressure 

(2,18)         (a^a)*  =aUa2+^b^+|  -^ 


Substituting  (2.16),  (2.17)  in  (2.18),  we  obtain 

6»  The  minimum  radius. 

The  value  of  a  quantity  q  at  the  time  when  the  bubble 
is  exactly  at  its  minimum  size  will  be  denoted  by  q.   The 
mlnlinum  radius  a  of  the  bubble  is  obtained  by  setting 
a  =  0   in  (2.17) •   The  resulting  equation  is  simplified  by  intro- 
ducing, in  place  of  a,  the  (non-dimensional)  internal  energy 
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u 

of 

the 

gas. 

defined 

by 

(2, 

.20) 

U    : 

k 

» 

•' = =(1)*" 


The  quantity  u  represents  the  portion  of  the  total  energy 
of  the  system  which  exists  as  internal  energy  of  the  gas  at 
the  time  of  minimuin  size.   Proa  (2.17),  we  get  the  following 
equation  which  determines  u: 

—      —  2 

(2.21)  -^^  =  I  {^    ,  (Graphed  in  Fig.  5). 

The  values  of  various  other  quantities  at  the  time  of 
minimum  size  of  the  bubble  are  of  interest.  From  (2.16)  we 
have 

(2.22)  I     =l^  =  l|uS 


a"    k 


and  (2.19)  yields 


8/3 
(2.23)         (a^a)   =   g/g  u    (4  -  3u) 


3 
81? 


Thus,  given  s,  the  internal  energy  u  at  minimuin 
size  is  obtained  from  equation  (2o21).  All  the  other  quantities 
a,  t,  (a^a)'   can  then  be  found  from  (2.20),  (2.22),  (2.23). 

7.  The  pressure  pulse. 

Sections  5  and  6  give  a  cociplete  description  of  the 
motion  of  the  bubble  when  it  is  near  Its  minimtjm  size.  We 
shall  now  investigate  the  pressure  pulse  delivered  to  the 
surrounding  water. 
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The  potential  function  ^    for  the  flow  of  the 
surrounding  water  has  been  constructed  in  equation  (2.1). 
The  pressiire  in  the  water  can  be  obtained  from  Bernoulli's 
equation  as  follows: 

(2.24)  £  ^  1  (grad  4-  )'  -  (||  -  B'  If-)  =  0, 

where  P   is  the  excess  over  hydrostatic  pressure,  Z  =  R  cos  0, 
and  the  term  B'  ^%-   appears  because  of  the  moving  coordinate 
system.   A  substitution  of  (2.1)  in  (2.24)  yields 

2 

(2.25)  ~=  ■     t    '        +  terms  in  higher  powers  of  ^  • 

For  points  in  water  not  too  near  the  bubble,  the  first  term 
on  the  right  hand  side  of  (2»25)  is  dominant  and  the  other 
tei'ms  may  be  dropped.   Introducing  non-dimensional  variables, 
we  obtain 


(2.26)         p.I^u!ai:,2Pa„  L(,2.,.. 

C^   R      3    R 


The  pressxn^e  in  the  water,  therefore,  depends  essentially 
on  the  quantity  (a  a)   • 

o 

The  forinula  (2.18)  for  (a  S)*  can  be  given  an  inter- 
esting  interpretation.   The  last  term  —   ..  ^ /■   represents 
the  (non-dimensional)  pressure  of  the  gas  inside  the  bubble, 
and  (2.26)  then  shows  that  the  pressure  at  any  point  in  the 
water  is  composed  of  two  parts:  the  internal  pressure  of  the 
gas,  and  the  dynamic  pressure  due  to  the  motion  of  the  water. 
This  is  more  clearly  seen  if  equations  (2.26),  (2.18)  are  re- 
written in  terms  of  the  original  dimensional  variables.   The 


CONFIDENTIAL 


C0Iv]FID3NTIAL 
-21- 


reault  is 


(2.27)  P  =  I   f|pA»2  +  ^pB«2  +  p(A)j       . 

where     p(A)      Is    the  pressure   of    the   gas   inside   the  bubble 
when  its   radius    is     A. 

The  principle  of  stabilization  asserts  that  the  most 
important  terra  in  the  expression  (2.27 )  is  the  gas  pressure 
terra     p(A).      To  obtain   the  maxiravun  possible   effect,    one 

should  maximize     p(A)    even  at   the   expense   of   complete   elinina- 

12    12 

tion  of  the  dynamic  terms  •$  PA'   +  -J  pB'   •   This  is  not 

obvious  and  requires  a  mathematical  proof. 

Returning  now  to  the  expression  (2»26)  for  P,  we 
wish  to  determine  the  time  when  the  pressure  pulse  reaches 
a  peak  value.   It  is  to  be  expected  that  the  peak  pressure 

occurs  when  the  bubble  is  of  minimum  size.  Although  this 

2»  • 
is  not  obvious  from  the  expression  (2.18)  for  (a  a)  ,  it 

is  a  sirjple  consequence  of  the  equation  (2ol9).  By  differ- 
entiating (2«19)  with  respect  to  a  and  by  using  (2.17),  we 
have 

d  f„2V«_    1  .  15  3"   11    k  _  :2    93^  ,21   k 

J3         2  *  * 

This  shows  that  ^(a  a)   is  negative  and  therefore  attains 
its  maximum  at  the  smallest  possible  value  of  a.   Hence 
the  peak  pressiire  occurs  at  the  time  of  the  minimum  size  of 
the  bubble. 

The  value  of  the  peak  pressure  P   can  now  be  obtained 
by  substituting  (2.23)   in  (2.26);  the  resulting  expression 
is 

P   Ij 

(2.28)  ^  =  ^  •  -~W5   •  ^^^^(4  -  3u)  . 
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This  Indicates  the  dependence  of  the  peak  pressure  P  on 
the  Internal  energy  u  at  the  time  of  minimum  size,  and 
through  (2. 21) J  on  the  linear  momentum  s. 

The  quantity  u  '  (4  -  3u)   will  be  called  the 
"pressure  factor"  and  will  be  denoted  by  the  symbol  q: 

0/3 
(2,29)  q  =  u   (4  -  3u)  . 

It  depends  on  s  by  virtue  of  (2.21). 
8»  The  optimum  peak  pressure* 

V/e  are  no-^.v  in  a  position  to  find  that  linear  momentum 
s  which  produces  the  maximum  peak  pressure.   The  peak  pressure 
P,  by  (2. 28),  is  proportional  to  the  pressure  factor  q  in 
(2»29),  which  in  turn  depends  on  i  by  virtue  of  (2»21).   A 
graph  of   q  as  a  function  of  t-^  is  drawn  in  figure  6,  and 
demonstrates  the  very  significant  fact  that  q  Is  largest 
when  s   is  practically  0  (or  u  practically  !)•  Actually, 
by  (2»29),  q  is  a  maxiimxm  when  u  =  ^  ,  or  by  (2.21),  when 
s  B  .151  k^.   But  since   k  is  small,  this  differs  so  little 
from  s  =  0  that  it  can  be  neglected.   See  figure  6. 

Thus,  the  following  general  Principle  of  Stabilization 
has  been  demonstrated:  For  a  given  roass  of  explosive,  the 
optimum  peak  pressure  in  the  secondary  pulse  is  obtained  by 
keepirg  the  bubble  motionless  at  the  time  of  its  minimum  size. 

The  value  of  the  optimum  peak  pressure   ^opt.   ^^^  ^® 
obtained  from  (2.28)  by  setting  u  =  1,  with  the  result 

p    -.1     ^o^ 

^opt,-  R  •    -^:^   • 


Using  (2»7)  and  (2.10),  and  expressing  P^  In  atmospheres 
by  means  of  P  =  D  /33  atmospheres,  where   D   is  the 
distance  in  feet  of  the  center  of  the  explosion  from  a  point 
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33  ft.  above  sea  level,  we  have 

(2, 30)  P^p^^  =  176  Jlj^  atmospheres 

where  W   is  the  weight  of  T.  TJ.  T.  in  pounds  and  R   is 
the  distance  from  the  center  of  the  bubble  in  feet. 

This  may  be  compared  with  the  peak  pressure  of  the 
primary  shock  wave,  whose  value,  experimentally  obtained,  is' 

^shock  ~  1160  — w —  atmospheres. 


The  optimum  peak  pressure  of  the  secondary  pulse  is  approxi- 
mately 15  percent  of  the  shock  wave  peak  pressure,  but  the 
duration  is  much  longer, as  will  be  seen  in  the  next  section. 

9.  The  impulse. 

The  impulse  I  per  unit  area,  carried  by  the  secondary 
pulse,  is  the  time  integral  of  the  pressure. 


I  =  /  PdT  =  C  /  Pdt, 


where  the  limits  of  integration  are  the  times  when  P  =  ot 
By  (2.26), 


T  _   o  •  L 
^  -  ~3~   ^ 


J   (a2a)'dt. 


2*  • 
where  the  limits  of  integration  are  the  times  when  (a  a)  =0. 

Using  the  fact  that  the  motion  is  approximately  symmetric  about 


*  See  [3]  page  13, where  this  experimental  result  is  quoted 
in  metric  units. 

^. 

^{.  The  significance  of  this  value  of  the  impulse  will  be 

discussed  later. 
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the   time  of   the  rainimujA  size  of   the  bubble,    the  above 
Integral  becomes 

(2.51)  I  =  -^   •  ^  .  [a'^a]. 


>•   o 


where  [a  a]   is  to  be  evaluated  at  the  time  when  (a  a)  =  0, 
or   a  a   is  a  ma'ximuj.i« 

Equations  (2.8),  (2.9)  and  (2.16)  yield 

4*2       ,  1/4    3s^    4 
aa  =a-ka'   -  — n  -  a 

2a'^ 


and 


,  4*2.    T     k   ^  3s^    .3 
(a  a  )  =1 -TTT  +  -rr 4a   . 

4a"^/^   a*^ 


Setting  the  last  e:cpression  equal  to  zero  in  order  to  find 

2*  2 

the  maximum  of   (a  a)  ,  we  obtain  an  equation  for  a.   If 

both  k  and  s   were  equal  to  zero,  the  solution  of  this 
equation  would  be  a  =  .es.   In  a  more  realistic  case  when 

k  =  .2,  s  =  .06,  one  obtains  a  =  .61.   This  illustrates 

2» 
that  the  value  of  a  when  a  a   is  a  maximum  is  not  very 

sensitive  to  changes  in  the  momentum  s,  since   s   is  generally 

quite  small  in  actual  casesi  Using  a  =  .61  and  k  =  .2, 

the  value  of  [a  a]  is  .55. 

By  (2.7),  equation  (2.31)  becomes 

^o^''^     83   W  ' 
I  =  .73  — ^ —  =  *'Y  1^  — — —  atmosphere-secondso 

D  /     ^ 
o 

Taking  a  typical  case  of  a  depth  of  100  ft.,  so  that  D  =133, 
we  find  that 

(2.52)  I  =  .37  — - —  atmosphere-seconds o 


i'c     See  Part  Hi. 
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We  have  shovm  that  the  impulse  carried  by  the  secondary- 
pulse  is  practically  independent  of  the  linear  momentum  s 
of  the  bubble,  and  is  given  by  (2. 32). 

The  impulse  carried  by  the  primary  shock  wave  is  experi- 
mentally measured  to  be" 

^shock  ~  •-2  — S —  atmosphere-sacoiids. 

A  comparison  with  (2.32)  indicates  that  the  impulse  carried 
by  the  secoadary  pulse  is  three  times  the  impulse  carried  by 
the  shock  wave. 

This  conclusion  can  be  stated  in  an  approximate  manner 
in  terms  of  the  durations  of  the  pulses.   Using  the  relation 
between  the  peak  pressures  obtained  at  the  end  of  the  last 
section,  one  can  state  that  the  secondary  pulse  lasts  about 
eighteen  times  as  long  as  the  shock  wave. 

This  result,  and  formula  (2.32),  should  not  be  taken 
too  literally  but  merely  as  order  of  magnitude  statements. 
The  reason  for  this  is  that  the  impulse  carried  by  the  secondary 
pulse  was  computed  between  the  times  when  the  pressure  is  zero. 
This  includes  a  relatively  long  period  of  time  during  which 
the  pressure  is  low,  and  the  resulting  contribution  to  the 
Impulse  will  not  be  experimentally  noticed  and  will  not  have 
much  effect  in  damaging  structures. 


-"<•  [3],    p.  13 
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Part  III.  The  Mechanism  of  Stabilization  by 
Gravity  and  the  Sea  Bed« 

1,  The  exact  equations  of  motion. 

In  this  part,  the  effects  of  gravity  and  of  a  rigid 
wall  upon  the  motion  of  the  bubble  are  considered.   The 
linear  momentum  s   produced  by  these  factors  in  the  course 
of  one  complete  pulsation  of  the  bubble  Is  calctilated.   As 
was  shown  in  part  II,  it  is   s  which  determines  the  peak 
pressure  produced  by  the  bubble  at  the  time  of  its  minimum 
size.   An  equation  for  determining  the  best  location  of  the 
mine  can  be  obtained  by  setting  the  momentum  s   equal  to 
zero* 


Dft. 


33  ft 


sea    LeoeL 


y/////////y///////////y 
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Use  the  same  variables  A,B   as  in  part  II  and  let 
D  stand  for  the  (fixed)  distance  of  the  sea  "bed  from  a 
point  33  ft.  above  sea  level.   In  place  of  the  equations 
(2.2),  (2.3),  (2.4),  we  have""' 

(3.1)      2npA^  [(l+f^)A^  -  2f^x\B  +  (i+f2)B^J 


+  -TTA"^pg(D-B)  +  G(A)  =  E, 
o 

where   f  ,f-,fp  are  functions  of  the  ratio  A/B  which 

represent  the  influence  of  the  rigid  wall  on  the  motion  of 

4    3 

the  bubble,  — ir  A  pg(D-B)  represents  the  gravitational 

potential  energy  due  to  the  lack  of  water  in  the  space 
occupied  by  the  bubble,  and  g  is  the  acceleration  due  to 
gravity.   Explicit  expressions  for  f  ,  f,,  fg  are  derived 
in  [5],  and  will  be  reproduced  in  appendix  II  of  this  reports 

Introduce  non-dimensional  quantities  as  in  equations 
(2.5),  (2.6)  of  part  II,  sections  3  and  4.   The  energy 
equation  becomes 

(3.2)  a^Cd+f^)!^  -  2f3_ab  +  {l  +  f^)h^] 


B 
where   d  =  y  ,  b  =  -2-  and  B^   Is  the  initial  distance  of 

the  bubble  from  the  sea  bed. 

The  motion  of  the  bubble  is  determined  from  (3.2) 
and  one  of  the  Lagrangean  equations  (2.12),  (2.13).   Using 


•3:-  See  appendix  II  for  details. 
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tJae  b-oquation    (2.13),   we   obtain 
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3,1 


3^   ! 


(3.3)  ^  [2a'^(^+f2)b     -  2a^f-^a] 


a 
2? 


df 


do6 


if. 


dT 


-  2  — =  ab   + 

d06  dOL 


2^2 


^^. 


where 
(3.4) 


oc   = 


a 


^   ' 


and  the   quantities     f 


are  fvuictions    of     06     only. 


The   differential   equations    (3.2),    (3.3)    are   to  be 
integrated,    subject   to   the   following   initial  condition: 


at      t  =  0, 


a  =  a. 


0 


b  =  b^,   b  =  0, 

4/3 
where   a   Is  the  smallest  root  (near  k  -^  )  of 
o 

a*^  +  —^jT  =   !•   Of  course,  these  Initial  conditions  are  not 

exactly  realistic,  since  at  the  very  beginning  a  shock  wave 
is  formed  by  the  explosion.   But  the  time  interval  required 
for  incompressible  flow  to  set  in  is  relatively  minute  and 
may  be  ignored. " 

The  quantity  in  the  brackets  on  the  left-hand  side 
of  (3.3)  is  the  linear  momentum  s   of  the  system.   The 
first  term  on  the  right- hand  side  of  (3.3)  is  due  to  the 
presence  of  the  rigid  wall,  and  the  second  term  is  due  to 

gravity. 

By  integrating   equation   (3,3)   v/e   obtain   the  ruomentvuti 
s      at   the   end  of   the  period  of  pulsation   of   the  bubble: 


*  See   [4], 
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(3.5) 


Jo     2b  "^ 


^o    *2        ^   ^1   •;         ^2    -2 


dt 


o      i/o 


3^4. 

a   dt 


where  the  integration  is  extended  over  the  full  period 
t   of  the  pulsation, from  the  explosion  to  the  time  of 
minlrfium  size» 

2«  The  approximate  evaluation  of  the  period  and  the  momentum* 

The  results  of  a  numerical  integration  of  equations 
(3«2),  (3«3)  are  tabulated  for  special  cases  in  appendix  II. 
They  serve  as  a  check  on  the  approximations  which  will  now 
be  made  to  evaluate  the  period  t  and  the  momentum  s» 

The  bubble  expands  to  a  maximiin  size  before  contracting 
again.   Indicate  the  value  of  a  quantity  at  the  time  of  maxi- 
mum size  by  a  subscript  1.   Thus   t-,   is  the  time  of  maximum 
size,  s-   is  the  linear  momentum,  etc.   We  shall  introduce 
the  follov;ing  approximations  which  are  especially  accurate 
when  the  bubble  is  in  its  balanced  position: 

1.  The  time   ^  and  linear  momentum  s   at  the  minimum 
size  of  the  bubble  is  twice  the  corresponding  quantity  at 

the  maximum  size;  i.  e.,  t  =  2t-,  s  =  2s,.   This  assumption 
agrees  very  closely  v/ith  the  numerical  integration  of  the 
equations.   It  means  that  the  motion  is  approximately 
symmetric  about  the  time  t^   of  maximum  size.   This  would 
be  exactly  correct  if  the  b-coordinate  did  not  vary. 

2.  During  the  first  half  of  the  period  of  pulsation, 
until  time   t-,,  b  is  small  and  can  be  taken  as  zero,  so  that 
b  remains  equal  to  b  .   This  agrees  satisfactorily  with 
the  numerical  integration,  as  well  as  with  experimental 
evidence. 
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Substltuting  b  =  "b^  In  (3.2),  we  get 


(3.6) 


e 

a  = 


The  solution  of  (3.6)  is  immediately  obtained  by  integra- 
tion, with  the  following  result  for  the  half -period  t^: 


(3.7) 


da 


where  a   is  the  smallest  and  a,   the  largest  root  of 
o  X 


(3o8) 


1  -  a*"  -  k  a 


-V4  =   0. 


To  evaluate  the  integral  (3.7)  it  is  necessary  to 
have  a  specific  value  for   k.  The  value  of   k  is  given 
by  (2.10)  and  depends  on  the  distance   D   of  the  explosion 
from  a  point  33  ft.  above  sea  level.   The  following  table 
is  calculated  from  (2.10): 


^0 

33 

83 

133 

183 

k 

.146 

.183 

.206 

.223 

Selecting  a  depth  of  water  of  about  100  ft.  as  typical, 
the  value  of   k  is  approximately  .2.   We  shall  select  this 
value  of   k  throughout  the  remaining  calculations.   A  differ- 
ent value  for   k  will  change  the  formulas  slightly. 

For  k  =  .2,  the  roots  a  ,  a,   of  (3.8)  are 


(3.9) 


a^  =  .118, 


a^  =  .924. 
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The   last   equation    (3.9)    is   noteworthy,    since      a^      repre- 
sents   the   non-dimensional  maximum  radius    of   the  bubble. 
The  actual  maximum  radius     A,      is   thus 


A^  =  .924   L. 

k 
If  the  Internal  energy  ■  ■„ /j   of  the  gas  were  neglected, 

the  maximum  radius  would  be  exactly  L« 


The  quantity  ^1  +  f^  in  the  integral  (3«7)  is 

very  closely  equal  to  1  +  ^=1+^   as  is  shown  in 

2       4bo 

appendix  II.   Making  this  approximation,  we  can  evaluate 
the  Integral  (3.7)  numerically.   This  numerical  evaluation 
Is  discussed  in  appendix  III,  the  result  being 


t^   =  .735  +  ^^ 


The  total  (non-dimensional)  period  t  of  the  pulsation  is 
therefore 

(3,10)     t  =  2t3^  =  1.47  +  i|^  =  1.47  I  1  +  i^  j   . 


Formula  (3.10)  shows  that  if  no  rigid  walls  are 
present   (b  =  oo),  the  actual  period  T  in  seconds  would 
be 

T  =  1.47  C, 

or  about  3/2  times  c  •   This  provides  the  physical  inter- 
pretation of  the  scaling  factor  C  mentioned  in  section  3 
of  part  II. 

Substituting  b  =  b^^,  b  =0  in  (3.5),  we  obtain  for 
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the  half -momentum  s^. 


or 


^  „        , df  ^ 

(3.11)  St  = ^   /    ^ da 


^     2b2  t/a„ 
o    o 


Vl  +  fo 


,  /--.  -^^^^^  ^. 


/»a^    a"/' 
%  ir. 


3    ,  -3/4 
a  -  ka  ' 


using  (3o6)«   The  first  term  on  the  right  hand  side  of 
(3.11)  is  the  downward  momentum  due  to  the  rigid  wall,  and 
the  second  term  is  the  upward  momentum  due  to  gravity* 
Again  taking  k  =  .2,  using  the  approximation 

; a    df  /dot  a 

Jl  +  f  =  1  +  ^  ,   ^'  =  1  -  5^  ,  and  evaluating 

these  integrals  numerically,  we  get  the  following  expression 
for  the  full  momentum  s: 

(3.12)  ^--  (^    -^)-  d:^  l-'''*^7r 

See  appendix  III  for  the  numerical  evaluations. 
3«  The  stabilized  position. 

By  the  principle  of  stabilization,  the  maximiAm  peak 
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pressure  is  obtained  by  setting  s  =  0,  with  the  result 


.704      + 

.148 
bo 

"°           .113 

b2 

0 

•  019 

Expressing  d   in  descending  powers  of  b  ,  we  get 


(3.13)         d  =  6.2  b^  +  3.3  b   +  .4   . 


A  graph  of  equation  (3.13)  is  drawn  in  figure  2. 

More  generally,  the  peak  pressure,  as  a  function  of 
the  distance  b   of  the  mine  from  the  sea  bottom,  is 
obtained  by  combining  equation  (2.2B),  (2.29),  figure  6, 
and  equation  (3.12).   As  a  typical  example  of  this  dependence, 
consider  the  case  discussed  in  section  3,  part  I,  of  a  mine 
containing  1500  lbs.  of  T.  N.  T.  in  water  of  depth  150  ft. 
The  scaling  factor  L,  as  obtained  from  figure  1,  is 

L  =  27.6  ft. 

The  depth  D  of  the  sea  bed  from  a  point  33  ft.  above  sea 
level  is  183  ft.,  so  that 

d  =  ^  =  6.63. 

The  best  location  of  the  mine  according  to  figure  2  is 
b^j  =  .76,  or 

B^  =  Lb„  =  21  ft. 
o     o 

above  the  sea  bed. 
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If  the  mine  is  located  at  various  other  distances 
B  ,  the  pressure  factor   q   can  be  obtained  from  figure  6  by 
first  determining   s   from  equation  {5,l'^)^      The  results  are 
tabulated  below,  and  the  graph  is  drawn  in  figure  3, 


^0 

.6 

.7 

.76 

.8 

«9 

1.0 

1.1 

1.2 

s 

-.066 

-.020 

+.003 

.015 

.039 

o059 

.074 

.087 

q 

.51 

.92 

1.0 

.96 

.72 

.55 

•  47 

.42 

As  stated  in  section  3,  part  I,  figure  3  shows  a  remark- 
able sharpness  in  the  peak  pressure  curve  as  a  function  of  the 
distance  from  the  bottom. 

4.  The  migration  of  the  bubble. 

A  forintula  for  the  distance  travelled  by  the  bubble  in 
the  couj?se  of  its  pulsation  is  difficult  to  obtain  because  it 
involves  a  complicated  repeated  integration.  But  by  combining 
a  theoretical  argument  with  the  results  of  numerical  integra- 
tion, an  empirical  f orimila  can  be  developed* 

It  seems  reasonable  to  suppose  that  the  displacement 

A  b  =  b  -  b   of  the  bubble  is  an  odd  function  of  s.   It 
o 

can  therefore  be  represented  by  the  beginning  terms  of  a  Taylor 
expansionj 


Ab  =  c,s  +  CgS 


where  c-,,  Cp  are  appropriate  constants.   In  fact,  a  theoretical 
Justification  of  this  can  be  given  on  the  basis  of  the  differ- 
ential equations  (3.2),  (3.3),  but  this  will  be  omitted  here. 


We  can  determine  the  constants   c 


empirically  by 


using  the  results  of  the  numerical  Integration  tabulated  in 
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appendix   I.      The    least   squares   solution  yields    the  following 
formula : 


(3.14)  Ab   =   19s    (1   -   62s^). 


A  graph  of  (3.14)  is  drawn  in  figure  7,  and  is  compared  with 
the  tabulated  values.  The  formula  (3.14)  is  seen  to  be  very- 
close  to  the  tabulations. 

5.  The  correction  due  to  the  free  surface* 

For  the  sake  of  completeness,  we  shall  include  some 
remarks  concerning  the  effect  of  the  free  water  surface.   The 
latter  exerts  a  repellent  force  on  the  bubble,  aiding  the 
downward  force  of  the  sea  bed.   A  complete  discussion  would 
require  a  knowledge  of  the  potential  function   <|>  and  of  the 
corresponding  values  for   f   f^  and  fg.   A  careful  examina- 
tion of  the  previous  argument,  however,  shows  that  only  the 

first  terms  in  the  expansions  of  f  and    o  in  powers  of 

°      ab 
a  were  used.   In  appendix  II  these  first  terms  are  obtained  by 

successive  reflections.   The  results  are  as  follows: 

If  C  is  the  depth  of  the  bubble  below  the  free  svirface, 

and  X   is  set  equal  to  z^^  ,  then  the  expansion  of  f   begins 


C+B 


with 


a 


^^'-' 


3f 


o 


and  the  expansion  of  -^^  begins  wi 


th 


-     -AyG(x) 
Sb'^ 

where   the   formulas   for  P(x)    and  G(x)    are   given  in  appendix   II. 
A  short  table   of  their  values  follows: 
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X 

1 

.9 

.8 

.7 

.6 

.5 

.4 

.3 

.2 

.1 

0 

P(x) 

1 

.86 

•  72 

.58 

.43 

.28 

.11 

-.06 

-.25 

-.46 

-.6  9 

G(x) 

1 

1.00 

1.00 

1.01 

1.03 

1.07 

lolO 

1.22 

1.33 

1.54 

1.86 

<33  ft.  aboue  sea-leoel 


sea  -  ieoel 


D 


B 


777777/F77777777/7777, 


7/777 


When  these  expansions  are  used,  forinulas  (3.10)  and 
(3. 12)  for   t  and  s   are  modified,  as  follows: 


(3.15) 


t  =  1.47   1  + 


.185  F(x) 


(3.16) 


s  =  - 


G(x) 


.113  - 


.019  P(x) 


*^ 


.704  +  '^^\^^^^ 
^o 
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It  is  to  be  noted  that  x  =  1  corresponds  to  no 
free  svirface   (c  =  ^'o  ),  while  x  =  0      corresponds  to  a 
point  midv/ay  between  the  free  surface  and  the  rigid  bottom. 

An  interesting  result  is  that  F{x)  =  0  when  x  =  l/3j 
so  that  the  influence  of  the  free  surface  on  the  period  of 
the  bubble  cancels  the  influence  of  the  rigid  bottom. 
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Appendlx  !•  The  Numerical  Integration  of  the 
Differential  Equations* 

To  test  the  validity  of  the  approximations  made  in 
part  III,  the  exact  differential  equations  (3.2),  (3.3) 
were  Integrated  nxmerlcally.   The  coefficients   f^,  f^,  fg 
and  their  derivatives  with  respect  to  ex  ,  are  tabulated 
in  table  4,  appendix  II.   The  ntunerlcal  integration  was 
carried  out  by  the  Mathematical  Tables  Project  operating 
under  the  Applied  Mathematics  Panel." 

In  this  report  we  have  not  reproduced  the  com.plete 
tables,  but  have  Included  the  graphs  drawn  in  figures  8  -  11, 
which  are  based  on  these  tables.  Here  we  shall  set  down  only 
the  most  Interesting  items  in  these  calculations,  namely,  the 
behavior  at  the  beginning,  at  the  time  of  maximum  size,  and 
at  the  time  of  minimum  size.""  The  following  table  shoves  the 
results: 


•55-  In  particular.  Dr.  G.  Blanche,  Dr.  C«  Lanczos  and 

Dr.  A.  N.  Lowan  helped  overcome  the  considerable  difficulties 

in  the  numerical  integration. 


-;:-5:-  Details  of  the  computation  can  be  obtained  from  an  extensive 
report  by  the  Mathematical  Tables  Project  prepared  for  the 
Applied  Itothematlcs  Panel. 
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Table   2,    The   results   of  nmnerlcal   Intei^ratlon 


t 

a  ' 

• 
a 

b 

• 

b 

s 

s 

I.    d  = 

11.17 
1.5 

0 
.8289 

.1172 
.925 

0 
0 

1.5 
1.565 

0 
.066 

0 
.0189 

.00017 
.079 

1.657 

.192 

0 

2.11 

16.0 

.0378 

.0006 

2.    d  = 
b     = 

0 

11.17 
1.17 

0 
.8505 

.1172 
.926 

0 
0 

1.17 
1.243 

0 
.0231 

0 
.00716 

.00016 
.0794 

1.6993 

.138 

0 

1.49 

15.8 

.0139 

.0 

3.    d  = 

11.17 
1.08 

0 
.8576 

.1172 
.926 

0 
0 

1.08 
1.156 

0 

.006 

0 
.00205 

0 
.0791 

1.7146 

.120 

0 

1.19 

6.9 

. 00393 

0 

4.    d  = 

11.17 
.9 

0 
.8753 

.1172 
.926 

0 
0 

.9 

.986 

0 

-.041 

0 

-.0154 

0 
.0769 

1.7574 

.167 

0 

.352 

-16.9  ' 

-.0270 

-.076 

5.    d  = 

2.0 

0 
.8043 

.117 
.922 

0 
0 

2.0 
2.01 

0 
-.0470 

0 
-.0127 

0 
-.0001 

1.610 

.163 

0 

1.58 

-17.4 

-.0251 

-.00236 

6.    d  = 

oo 

0 

.117 

0 

1.5 

0 

0 

0 

^= 

1.5 

.8204 

.923 

0 

1.52 

-      .0769 

-.0219 

-.0002 

1.6504 

.202 

0 

.808 

-15.4 

-.0426 

-.007 

7.    d  = 

oo 
1.2 

0 
.8266 

.117 
.922 

0 
0 

1.2 
1.23 

0 
-     .111 

0 
-.0339 

0 
-.0004 

1.699 

.237 

0 

.282 

-12.6 

-.0645 

-.58 

•M-  Tbls    case   corresponds    to  a  rigid  wall  alone,   with  gravity  ignored. 
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The  period  t  and  the  momentum  s   for  the  cases 
listed  above  were  computed  from  the  formulas  developed  in 
part  III.   They  are  collected  in  the  following  table,  and 
compared  with  the  numerical  values  given  in  the  preceding 
table*   The  agreement  is  excellent. 


TABLE  3, 


t 

■ 

c 

- 

Calculated 
from  (3.10) 

Numerical 
(Table   2) 

Calculated 
from   (3.12; 

Numerical 
(Table  2) 

1.    d  = 

11. 17 , 

b^=  1.5 

1.55 

1.66 

+•039 

+.038 

2.    d  = 

11.17, 

b^=  1.17 

1.70 

1.70 

+.012 

+.014 

3.    d  = 

11.17, 

b^=  1.08 

1.72 

1.71 

+.002 

+.004 

4.    d  = 

11.17, 

V    -9 

1.77 

1.76 

-.030 

-.027 

5.    d  = 

oo    , 

b^=  2.0 

1.61 

1«61 

-.026 

-.025 

6.    d  = 

oo    , 

^0=  1-^ 

1.65 

1.65 

-.044 

-.043 

7.    d  = 

oo    , 

b^=  1.2 

1.70 

1.70 

-.067 

-.065 
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Appendix  II.  The  Velocity  Potential  for  a  Pulsating 
Sphere  Moving  Perpendicularly  to  a  Wall> 

!•  Statement  of  the  problen* 

An  ideal  inconipressrole  fluid  is  bounded  by  a  plane 
infinite  rigid  wall^'  and  has  a  pulsating  and  moving  sphere 
imp.iersed  in  it.   It  is  required  to  find  the  velocity  poten- 
tial describing  the  (Irrotational)  flow,  and  to  obtain  the 
kinetic  energy  of  the  liquid.   The  case  of  a  moving  sphere 
is  classical  and  can  be  found  in  the  standard  books  on  Hydro- 
dynamics.  The  pulsating  sphere,  however,  does  not  seem  to 
be  generally  known  and  has  certain  features  of  interest. 
The  treatment  here  follows  [5]« 


Let  A  be  the  radius  of  the  sphere  at  any  time  T, 
and  B  the  distance  of  its  center  from  the  rigid  wall.   The 
velocity  potential  ^l*  describing  the  flow  of  the  liquid 
satisfies  the  potential  equation 


^-^^   ^xx    yy    zz 


«  The  method  used  here  can  also  be  applied  to  a  spherical 
wall* 
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and  has  the  following  boundary  conditions: 
(1) 


^   -  A'  -  3'  cos  0   on  the  sphere 
o  n 


-  ^-3-  =  0  on  the  plane, 

3  n  r-     7 


where   n  is  the  normal  direction  pointing  into  the  liquid 
and  the  "prime"  denotes  a  time  derivative* 

The  construction  of   <$>  can  be  decomposed  into  two 
simpler  problems  by  setting 


(2)  4>  =  A'  «|>^  -  B'  i^ 


where  <$»  ,  ^ -,      are  potential  fxinctions  satisfying  the 
following  boundary  conditions: 


(3)   -    Q  =  1      on  the  sphere,  -  ^ — -^  =0  on  the  plane 


3  $  X  3  ^ 

(4)   -  =  cos  0  on  the  s-ohere, ^ — 1  =  0  on  the  plane. 

3  n.  "         <^  ^ 


Physically,  <|>   represents  the  potential  function  for  a 
sphere  of  fixed  center  and  radius  expanding  with  unit  velocity, 
while  ^^      is  the  function  associated  with  a  rigid  sphere 
moving  away  from  the  wall  at  unit  velocity.   The  functions 
^  ,  ^ 1   will  be  constructed  by  means  of  the  method  of  images. 

2»  Some  theorems  on  im^xges* 

The  theorems  that  follow  supplement  well  known  theorems 
on  images.    See  Milne -Thorns on.  Theoretical  Hydrodynamics, 
Chapter  15« 
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We  wish  to  find  a  potential  function  <^  ,  defined 
outside  a  f-iven  sphere,  with  a  given  distribution  of  sources 
outside  the  sphere  and  vanishing  normal  derivative,  -^-   =0, 
on  the  "boundarjo   A  standard  method  for  constructing   <|>  is 
to  place  a  suitable  arrangement  of  sources  inside  the  sphere, 
called  the  image  of  the  original  distribution. 

The  image  of  a  point  soujrce  and  a  radial  dipole  are 
well  known,  and  will  be  found  In  Mi Ine -Thorns on  ,  pp.  420,421. 
Here  we  consider  the  image  of  a  line  so\ircec 

Theorem  1«  Consider  a  sphere  with  center  0  and  radius 
a.   The  image  of  a  uniform  line  source  of  strength  ix,     per  unit 

length  stretching  from  Q-j^   to  Qg   Is  the  following:  a  uniform 

OQp 
line  source  of  strength  M— g—  per  unit  length,  stretching 

between  the  points  Qn,Qo  inverse  to  Q-, 

line  sink  of  strength 


Qn,Q2   inverse  to 
QTQo 


Qrj,  and  a  uniform 


ar 


1'  "^2' 
per  unit  length  extending  from 


0  to  Q^« 
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Proof.   If  P   is  any  point,  the  Stokes  stream 
function  "^  (see  Milne-Thomson,   Chapter  15,  esp.  pp.  420, 
421)  is 

OTQp      ,         ^2  ,  '   'v 


'      a 
If  P   is  on  the  sphere,  we  have  r^  =  r,  , 

r„  =  — ^  r„  (  a  property  of  inverse  points,  due  to  the 

similarity  of  triangles   OPQ^  ,  OQ^P),  and  consequently 

-\j/'  =  -  ix  — i-^  =  const.   Therefore    -^-   =0  on  the  sphereo 
Combining  this  result  with  the  known  image  of  a  point 
source  (Milne -Thorns on,  pp.  420,421),  we  obtain  the  following; 

Theorem  2.  The  im.age  of  a  point  source  of  strength  m 
situated  at  the  point   Qg,  together  with  a  uniform  line  sink 
of  strength  ^    per  unit  length  extending  from  Q,   to  Q^, 
where  m  =  MQ-iQp   (so  that  the  total  strength  is  zero),  is 

the  following:  a  point  source  of  strength  m  -— —  at  the 

point  Qp   and  a  uniform  line  sink  of  strength   A*- "~a   ^^^ 

I    t 

unit   length  extending  betv/een     Q,,    Qg* 

5*    The   construction   of    <|>o» 


c 


The  function   ^    is  required  to  satisfy  the  boundary 
onditions  (3).   If  the  plane  were  not  present,  ^        would  be 


o 


a2  2 

^  ,  which  represents  a  source  of  strength  A   placed  at  0. 

The  presence  of  the  plane,  however,  causes  the  boundary  condi- 
tion (3)  on  the  plane  to  be  violated.   To  satisfy  it,  introduce 
the  image  of  the  source  relative  to  the  plane,  which  is  a 
source  of  equal  strength  at  the  reflected  point  Q^.   But  the 
boundary  condition  on  the  sphere  is  nov/  violated,  and  to 
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remedy  it.    Introduce   the    image   of    the  point   source   at     Q- 

relative  to   the   sphere.      This    image    is    a  point   source   of 

2  A  ' 

strength  A  ■=—  at  the  inverse  point  Q,   and  a  uniform 

line  sink  of  strength  A  x  P^r  unit  length  extending  from 
0   to   Q, ;  this  process  continues  and  requires  the  use  of 
theorem  2  for  the  images  relative  to  the  sphere*  The  successive 
Images  fit  very  neatly  together,  as  indicated  in  the  diagram. 
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The  successive  reflections  are  described  in  the 
following  table: 


Strength  of 
point  source 

Situation  of 
point  source 

Strength  per  unit 

length  of  uniform 

line  sink 

Situation 
of  uniform 
line  sink 

a2 

0 

a2 

^1 

A^  A 
0^1 

t 

A 

oq; 

a2  " 

^2  A   A 
0«1  0«2 

«2 

T 

A 

^  A 

etc« 


Set 


(5) 


and 


(6) 


^n  =  ^1-^2 ^n 


0^ 
By  the  inversion,  d   is  also  equal  to  — x-  •   The  general 

expressions  for  the  strengths  of  the  sources  are: 


(7) 


A  point   source   of  strength     AD       at      Q^     and     Q„.t? 
a  vuaiform  line   sink  of  strength     =r-     per  unit   length 
along     q\^^^     and     VlVs' 
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Introduce  the  quantity 


<«'  '^'^ 


The   quantities      d       can  be   expressed  in    terras   of    C5C  • 
We  have 


,  _        A        _  A  _  1 


^■^^      ^  2R  -  "no'       ^^       T         * 


Thus 

(9)  d^     =   oc,...,   d^^^  =  Y 


06 


"^1 


SO  that  the   d's   are  convergents  to  a  continued  fraction. 

The  quantities   D   can  be  obtained  as  functions  of 
Ci^  from  the  definition  (6).   They  can  also  be  obtained  by 
solving  a  second  order  linear  difference  equation  with 
constant  coefficients  for  the  quantity  =r-  ,  obtained  from 
the  rec\arslon  formula  (9)»  ^ 

It  is  convenient  to  expand  <|    In  the  neighborhood 
of  the  surface  of  the  sphere  in  terms  of  spherical  harmonics. 
The  expansion  of  a  unit  point  source  lying  outside  the  sphere 
at  a  distance  X  from  the  center  is 

2 
(10)        ^  "^  ^  P]^(cos  G)   +  ^  P2(cos  e)  +  ...   , 

X  X 


while  a  point  source  lying  inside  the  sphere  has  an  expansion 

2 
(11)        ^  +  -^P^  (cos  0)   +  ^  P2(cos  Q-)    +  ...   , 
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where  {p  (cos  9)|  are  the  Legendre  polynomials  of  cos  9© 
(See  diagram  on  page4i'.   Integration  with  respect  to  X 
yields  the  expansion  of  a  uniform  line  source.   The  ex- 
pansion of  "t    is  then  easily  obtained  by  using  the  distri- 
bution of  sources  and  sinks  given  in  (?)• 

For  the  present  purpose  it  suffices  to  know  merely 

the  value  of  — =H?  1^   dS,  the  mean  value  of  ^^  over  the 

4TrA2^  °  ° 

surface  K  of  the  sphere.   The  integral  of  a  unit  point 

so\xrce  over  the  surface  of  the  sphere  is,  by  (10)  and  (11), 
(12)  ^    or    i  , 

according  as  the  source  lies  outside  or  inside  the  sphere. 
The  first  of  these  is  to  be  expected  from  the  mean  value 
theorem  for  potential  functions,  while  the  second  is  a 
constant,  independent  of  the  position  of  the  sources  inside 
the  sphere. 

The  value  of  — ^^  f  ^^^      can  now  be  obtained  by 

using  the  distribution  of  sources  and  line  sinks  (7).   The 

contribution  of  the  distribution  ins  Ida  the  sphere  reduces 

PI  2 

to  A  -^  due  to  the  source  of  strength  A   at  the  center. 

This  follows  from  the  second  result  in  (12)  and  the  fact 

that  the  successive  images,  each  consisting  of  a  soiirce  and 

a  line  sink,  have  a  total  strength  zero.  The  contribution 

due  to  the  source  AD   at   0^+1  ^^'^   ^^^  line  sink  of 

strength       per  unit  length  along  ^^4.1  i^^^   equation 

n-1 
(7)    )    is 

«2.        1  A      ,        ^ 

n+1 

or 

«    T>                 At  n+1 

A  D„  .  n    -  tp: log 


"■^^       ^n-1        "     ^i 
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Thus   v/e   obtain  the   final  result 


(13) 


h 


dS  =  4TrA' 


1    +    Dn 


oo 

n  =  i   \ 


The  term  D,   is  due  to  the  point  source 
4,  The  construction  of  <!►•,• 


n+1   D. 
2 


A 


n-1 
at   Q 


log 


x) 


1* 


The  function  <|>  ,   is  required  to  satisfy  the  boundary 

conditions  (4),   If  the  plane  were  not  present,  <|>  ^   would 

be  the  velocity  potential  due  to  a  dipole  of  strength  ^ 

A   cos  0 
placed  at  0,  namely,  -« -^ —  •   The  boundary  condition  on  the 

plane,  however,  is  violated  and  to  satisfy  it,  introduce  the 

A  3 

image  radial  dipole  of  strength  -  %-  situated  at  the  reflected 
point  Qti  (see  the  diagram  on  page  45  )•  To  remedy  this  boundary 
condition  on  the  sphere,  which  has  now  been  violated,  introduce 
the  image  of  this  radial  dipole  relative  to  the  sphere*  This 
image  is  another  radial  dipole  of  strength 


3^ 

inverse  point  Q'.  (See  Milne-Thomson,p.  421  ).Etc« 


at  the 


The  successive  reflections  are  described  in  the  following 


table; 


Strength  of 
Dipole 


A^ 
2 

•i 
2 


A^     A^ 


a3_^ 


A^     A^     A^ 


Situation   of 
Dipole 


0 


etc* 
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The  general  expressions  for  the  strengths  of  the  dlpoles 
are: 

at  (3   ,  there  is  a  radial  dipole  of  strength  -^  Dv,  » 


(14) 


at  Qj,.n       II     II    tl     II         II      II         "  «   -^  T)^ 


3 

T   ^n  • 


The  expansion  of  a  unit  radial  dipole  lying  outside 
a  sphere  at  a  distance  X  from  the  center  is 

2 

(15)  -  A  -  2H  p^(cos  0)  -  §_  P2(cos  e) , 

while   a  radial  dipole   inside   the  sphere   has   an  expansion 

(16)  ip-,(cos   ©)    +^Pp(co8   ©)+ 

j^2    1  R"^     "^ 

These  are  obtained  from  (10),  (11)  by  differentiation  with 
respect  to  X»   See  Milne -Thorns on,  pp.  442,443. 

For  the  present  purpose  we  are  interested  in  the 
expansion  of  4>,   only  up  to  terms  involving  P,  (cos  9). 

Using  (14),  (15),  (16),  we  obtain 

oo      3  oo       3 

(l-^)    \  =  #  Z!  =4-     ^     ^^    Z!  =1-  R  Pi(cos  9) 


,3^3  P-,(cos  9) 

+  ~-  >  D  — = — 5 +  terms  in  higher 

nTo      R         Legendre  polynomials. 

It  is  convenient  to  evaluate  /  <i>  ^  dS  and  /  ^  ,  cos  9dS  . 
Using  the  orthogonality  property  of  Legendre  polynomials. 
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and   (17),    (5),    (6),   we  have 

(18)  /*  l^S   =  2nA=  Z  D^  4l 

v  K  in  =  ft 


'K  n  =  o 


(19)  f  ^  ^cos    9c3S   =|ttA^(1  +   ^^D^) 

5.    The   kinetic   energy  of   the   water* 

The  kinetic   energy    '^f    of   the  water   is 


dS 


by  Green's    theorem.      By   (2),    and  the   Green's   theorein 

d4> 


1/      o  an  J      1  a„ 


3n         ' 


7'=  -  i  /"[^'Ao^  ^  -  ^^/*i^  ^  *  '^fi^  ^ 


The  boundary  conditions  (3),  (4)  yield 

(20)  7^=  I/O U^/^o  ^  "  ^*'^/  *  i^S  +  B^/  $-L  CO 

L   ^K  '^K  I'K 


s  OdS 


The  integrals  appearing  in  (20) have  been  calculated 
in  the  preceding  sections.  By  (13),  (18),  (19), 

(21)*^=  2n^A^   (l+fo)A^  -  2f^AB  +  (^  +  f2)B^   , 
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where  the  coefficients  f  ,  f , ,  fp  are  the  following 

A 
functions  of  the  quantity  ^  =  ^  ! 


(22) 


oo 

n»  1 

Vi 

1      •,.-  *n+l' 

^i=IZ^^ 

4i 

^2=iZ"n   • 


n=  1 


The   expressions   for     d  ,    D       as   functions   of    ol   are 
given  in   (9)    and    (6), 

The   series   for     f   ,   f,,   fg      converge   for  oc   =  1/2 
(which  means     A  =  B)«      They  are  tabulated  in  table  4  below. 
Expansions   in  powers   of     oc    begin  as   follows: 


1      4 
f     =  oc  +  ^  ex.       +  •.. 

o  2 


^1  -  2*^      ''2°^     ■*■   ••• 


^2  =  2  "*■  2°^      ■*■   ••• 


V 


df./doc 


The   expansions   for       /T~+~f    ,        

°       /I  +  f( 
in  part   III,   begin  as  follows: 


required 


4 


l  +  f_^  =  H.|«. 


(24) 


af^d 


06 


/l^fo 


1-±C.+ 
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df 

dfn 

dfo 

oc 

^o 

^1 

^2 

0 

do. 

1 

da 

2 

da 

0 

0 

0 

0 

1.00 

0 

0 

.02 

.020 

.000 

.000 

1.00 

.02 

.00 

04 

040 

001 

000 

1.00 

04 

00 

06 

060 

002 

000 

1.00 

06 

01 

08 

080 

003 

000 

1.00 

08 

01 

.10 

.100 

.005 

.001 

1.00 

.10 

.02 

12 

120 

007 

001 

1.00 

12 

02 

14 

140 

010 

001 

1.01 

14 

03 

16 

160 

013 

002 

1.01 

16 

04 

18 

181 

016 

003 

1.01 

18 

05 

.20 

.201 

.020 

.004 

1.02 

.20 

.06 

22 

221 

024 

005 

1.03 

23 

07 

24 

242 

029 

007 

1.03 

25 

09 

26 

263 

035 

009 

1.04 

28 

11 

28 

283 

040 

Oil 

1.05 

30 

13 

o50 

.305 

.047 

.014 

1.08 

.33 

.15 

?2 

326 

053 

017 

1.09 

36 

17 

34 

348 

061 

021 

1.11 

39 

20 

36 

371 

070 

025 

1.15 

43 

23 

38 

394 

078 

030 

1.18 

48 

27 

.40 

.418 

.088 

.036 

1.23 

.53 

.32 

41 

430 

094 

039 

1.26 

57 

35 

42 

443 

100 

043 

1.29 

60 

38 

43 

456 

106 

047 

1.34 

65 

42 

44 

470 

112 

051 

1.40 

69 

46 

.45 

.484 

.120 

.056 

1.45 

.76 

.52 

46 

499 

128 

062 

1.54 

84 

60 

47 

515 

136 

068 

1.64 

95 

70 

48 

532 

147 

076 

1.81 

1.12 

86 

49 

551 

159 

085 

2.26 

1.55 

1.27 

.50 

.577 

.178 

.101 

CXD 

oo 

oo 

Table  4. 
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A  coirparlson  with  the  complete  table  shows  that 

1  +  -^  is  a  very  close  approximation  to  /l  +  f   for 

1 
all  oc  5  ■i,  and  differs  from  it  by  less  than  1  percent.   On 

the  other  hand,  1  -  -^  is  not  too  good  an  approximation 

(ii'f./doc                                  df^/dcx 
to  —      •     Further  terms  in  the  expansion  of 

/l"^  /I  +  fo 

are  not  necessary  since  they  produce  only  small  deviations 

from  the  forimlas  (3.12)  and  (3.11). 


6.  Approximate  theory  for  a  sphere  pulsating  between  a  rigid 
wall  and  a  free  surface. 

Consider  a  pulsating  bubble  at  a  distance  b  from 
a  rigid  wall  and  a  distance   c  from  a  free  surface.   For 
convenience,  the  free  surface  will  appear  vertical  in  the 
diagram. 


Q. 


Q. 


Q. 


to 


0 


% 


R 


R 


n 


If  0   is  the  original  position  of  the  center  of  the 
sphere,  the  successive  images  with  respect  to  the  rigid  wall 
and  free  surface  will  be  at  P,j  Q,;  Pp,  Qg;  P,,  Q,;  etc. 
Their  distances  from  0   will  be  given  by  the  following 


See  [6]. 
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foriflulas: 

J5T^  =  2m(b+c)  if     n  =  2ra 

=  2in(b+c)   +  2b      if     n  =  2m  +  1 

^     =  2m(b+c)  If     n  =  2iii 

=  2in(b+c)    +  2c      If     n  =  2in  +  1 

The  signs  of  the  images  are  given  in  the  diagramo 
Images  with  respect  to  the  sphere  will  be 
neglected. 

The  first  term  in  the  expansion  of  f   in  powers 
of   a   can  be  obtained  by  finding  the  potential  at  0   of 
the  Infinite  series  of  images  P.,  Q   Pp,  Q-,  etc.   Some 
simple  algebraic  transformations  lead  to  the  formula 


^o  =  HT^(2^^^^^  -  ^°e  2) 


where  x  =     i^  and 


1      1^1 


f(x)  =  p ^  +  p  ••• 

1-x^   9-x'^   25-x^ 


It  is  to  be  noted  that 

4r= ^  tf(^)  +  xf'(x)] 

Sh  (b+c)^ 

To  facilitate   comparison  with  previous   results. 
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f^     and     can  be  written  as   follov;s: 

°  Sb 


f ^  =  ^  (l-x)(2xf(x)  -  log  2)  =  I^F(x) 


3f  o 

— ^  =  -  -~  2(l-x)'^[f(x)  +  xf'(x)]  =  -  -^  G(x) 
o  b      2b^  2b'^ 


A  short  table  for  F(x)   and  G(x)   is  given  on  page  36^ 
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Appendljc   II I»   Nxunerical  Evaluation  of   Some 
Definite   Integrals 

m  part   III,    the  problem  of   evaluating   the    definite 
integrals    in  equations    (3.7),    (3,11)   arose.      These   integrals 
have   square  root   singularities   at   the   limits   of    Integration. 
In  this   appendix  vie   shall  develop  a   quadrature   formula,   which 
was  used  for  this  evaluation  but   does   not   seem  to  be   as  well- 
known  as   it   should.      It   is  based  on  Tchebycheff  polynomials. 

The   nth  Tchebycheff  polynomial     T   (x),   -i  ^  x  ^  1,  is 
defined  by 

T   (x)   =  cos(n  arc   cos  x). 


or  in  other  terras,  T  (x)  =  cos  nO  ,  where  x  =  cos  6. 

The  Tchebycheff  polynomials  are  orthogonal  with  respect  to 

the  weight  function  --rr=r  ,  i»  ©•, 

'  T  (x)  T  (x) 

^  ^   ■dx  =  0       ifn^m. 


^ 


^2 


This  follows  immediately  when  we  make  the  transformation 
X  =  cos  9. 

By  virtue  of  the  orthogonality  property,  the  follow- 
ing theorem  can  be  proved: 

Theorem.  Let  x^,  Xp,...,  x   be  the  n  zeros  of 
T  (x)  in  the  range   -1  ^  x  ^  1.   Let  f (x)  be  any  polynomial 
of  degree  at  most  2n-l.   Then 


,i/'-^ax=iZf(x,). 
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The  proof  is  similar  to  the  proof  of  the  Gauss 
quadrature  formula.   This  theorem  indicates  an  exact 
evaluation  of  the  integral  involving  2n  parameters 
(a  polynomial  of  degree   2n-l)  in  terms  of  n  specially 
selected  points.   It  is  exactly  analogous  to  the  more  well 
known  Gauss  quadrature  formula,  but  it  is  simpler  in  two 
ways:  the  zeros  x,,..o,x  of  T  (x)   are  easy  to  obtain, 
since  they  are  merely  the  zeros  of   cos  n  9  where  x  =  cos  9; 
and  the  weight  factors  multiplying  f(xi)   are  all  equal  to 
1/n. 

If  f(x)   is  not  exactly  equal  to  a  polynomial  of 
degree  2n-l  but  can  be  approximated  by  one,  we  can  write 
the  approximate  quadrature  fonmila 


f(x) 


'l(^- 


99 


2       ""  Til" 


The  accuracy  of  this  approximate  formula  depends  on  the 
closeness  with  which  f(x)   can  be  approximated  by  polynomials 
of  degree  2n-l. 

The  integrals  in  question  in  part  III  are  of  the  form 


a^ ^ 

—  da 


where  v  is  some  exponent,  and  a  ,  a_   are  the  two  zeros  of 
the  expression  under  the  square  root  sign.   This  can  be 
written  as 


I 


V   ^(a^  -  a)  (a  -  a^) 

•       ■  ■  '  "  '     •  — 
k 


da 

3 


1  -  ^"  -  -^   ^/(a3_-a)(a  -  a^) 
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\J{a     -  a)(a   -  a    ) 
where    the  function        a   •   —        ■  — —     has   no 

T  3  k 

singularities  and  Is  rather  smooth.   The  Integral  can  now 
be  evaluated  by  using  the  Tchebycheff  quadrature  forimila 
(changing  the  range  of  Integration  by  a  linear  transformation 
from  a  ,  a,   to  -1,  1).   The  integrals  were  evaluated 
in  this  way  by  using  the   5   zeros  of    Tg(x),  with  the  re- 
sults stated  in  part  III.   The  accuracy  depends  on  the  close- 

\/(a^  -a)  (a-a  ) 

ness  with  which  the  remaining  factor  a^. — 

/I    3    i^ 

can  be  approximated  by  a  9th  degree  polynomial.   It  is  fortunate 
that  the  Integrand  need  be  computed  only  for  5  Intermediate 
points.  (Also,  the  values  obtained  by  using  the  three  zeros 
of   T^(x)   agree  within  1  percent  with  the  values  quoted.) 

It  is  of  interest  to  see  how  the  momentum  s   depends 
on  the  parameter   k.   If  the  explosion. takes  place  near  the 
surface,  say  for  a  model  experiment,  the  value  of  k  would 
be  approximately  .16.   For   k  =  .16,   the  computations  yield 

in  place  of  (3.12).   All  the  constants  have  Increased,  but 
the  alternate  additions  and  substractlons  tend  to  cancel 
these  increases.   Thus,  the  new  equation  for  the  stabilized 
position,  where  s  =  0,  is 

d  =  6b^  +  3.32h  +  o4 
o       o 

which  is   practically   identical  with    (3.15).  The   total  period 
for   this   case  turns   out   to  be 


t   =   l.48/l   +  429| 
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wMch  is  practically  identical  with  (3.10). 

Finally,    if   the  internal  energy   is  neglected,    i.e., 
if     k  =  0,    the    Integrals  become  Beta  functions   and 

o 


I  =   1.49   1   +  1§23      „ 


The  stabilized  position  where   s  =  0  is 


d  =  5,lb^  +  3«lb„  +  04  . 
o       o 


These  deviate  somewhat  from  (3. 12)  and  (3»13)» 
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STOD/IARY  OF  FORMULAS 


!•      Scaling  factors    (for  T.    N.    T. ) : 


1/3 
L  =  13.2   |4rl  ft. 


C  =  2.85 


W 


1/3 


D 


5/6 


sec. 


where     W  =  weight   of  T.    N.    T.    In  lbs.,   D     =  distance.    In  feet, 
of   the  explosion  from  a  point  33  ft.    above   sea   level. 

2.      Period  of  pulsation: 

T  =  1.47C  =  4.19  "'    g,.;.    sec. 

D  ^6 


in  the  absence  of  outside  Influences;  If  the  bubble  Is  at  a 
distance  b   from  a  rigid  wall,  this  Is  to  be  multiplied  by 


1  +  1 


185 


3.  Maximum  radius; 

If  the  Internal  energy  of  the  gas  Is  neglected,  the 
maxlmuni  radius  of  the  bubble  is  L.   If  the  Internal  energy  is 

allowed  for,  the  maximum  radius  is  a,L,  where  a.   Is  the 

3    k 
largest  root  of  a  +  ■  i  j  =  1,  and  k  Is  given  by  the  formula 

k  =  .0607  D^^   . 

4.  Non-dimensional  momentum  at  time  of  minimum  size: 

The  momentum  due   to  gravity  and  a  rigid  bottom  is 
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-  1       /      TT-z  *019 


VAh'*"^] 


5.      Peak  pressvire  for   the    secondary  pulse: 
The  peak  pressure   is 


wV3 
176  — -r —  q  atmospheres 

R 


where  R   is  the  distance  in  feet  from  the  center  of  the 
bubble  and  q  is  a  factor  whose  graph  is  drawn  in  figure  6. 
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